Multicomponent fluid of hard spheres near a wall 
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The rational function approximation method, density functional theory, and NVT Monte Carlo 
simulation are used to obtain the density profiles of multicomponent hard-sphere mixtures near a 
planar hard wall. Binary mixtures with a size ratio 1:3 in which both components occupy a similar 
volume are specifically examined. The results indicate that the present version of density functional 
theory yields an excellent overall performance. A reasonably accurate behavior of the rational 
function approximation method is also observed, except in the vicinity of the first minimum, where 
it may even predict unphysical negative values. 
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I. INTRODUCTION 

Many interesting physical phenomena such as wet- 
ting and adsorption involve fluids and their mixtures 
at a solid-fluid interface. It is therefore not surprising 
that a lot of work has been devoted to these problems 
within the last three decades. The methods that have 
been used to study them include the integral equation 
approach of liquid state theory (see for instance Refs. 

jOToaches (see 
15|) and com- 

111 [13, El). 



1|, l^, 1^ 0, H, 1^,0, ISl ) ■ density functional a 
for instance Refs. H, i, [ip, [ul, [H, [H, [3, 

puter simulation (see for instance Refs. @, 

A rather simplifled but essentially correct physical pic- 
ture of adsorption may be obtained if one considers the 
solid surface as a planar smooth hard wall and the fluid 
as consisting of hard spheres (HS). Interestingly enough, 
this simple model is capable of accounting for the most 
important feature resulting from the interaction between 
the fluid particles and the wall, namely the strongly os- 
cillatory nature of the density distribution proflle of the 
fluid particles in the interfacial region. The particles are 
depleted from the surface of the wall due to excluded 
volume effects and this is the source of the density os- 
cillations. A particular realization of such a model is 
obtained from a binary HS mixture in which one of the 
species is taken to have an infinite diameter and to be 
in vanishing concentration. In this instance, the wall- 
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particle pair correlation functions lead immediately to 
the density profiles. In fact, the availability of the ana- 
lytical solution of the Percus-Yevick (PY) equation for 
additive HS mixtures obtained by Lebowitz allowed 
Henderson et al. [H to derive the density profile of a HS 
fluid near a hard wall within the PY approximation al- 
ready thirty years ago. More recently, Noworyta et al. 

also used the PY theory to study two binary mix- 
tures of HS near a hard wall. In the same paper they 
also performed grand canonical ensemble simulations for 
this system and used a version of density functional the- 
ory (DFT) to derive the density proflles of both species. 
Their results indicated that the second order PY theory 
provided in general the best agreement with the simula- 
tion results and that the DFT was a little less accurate. 
It is interesting to point out that for the most asymmetric 
case that they examined, that of size ratio 1:3, anoma- 
lies were observed, namely a huge discrepancy between 
theory and simulation, for which they could flnd no ex- 
planation. This work has served as a motivation for the 
present paper. Here we also tackle the problem of binary 
HS mixtures at a planar hard wall. In addition to the 
PY theory, we use the results of the rational function 
approximation (RFA) method ^2^, a different version of 
DFT 21], and NVT Monte Cario simulation. Apart from 
complementing the work of Ref . Q , our aim is to assess 
the value of both the RFA method and of our version 
of DFT to deal with this problem in a situation where 
both theories are subject to rather stringent conditions, 
namely when there is a disparate size of the diameters 
but both species occupy a similar volume. 

The paper is organized as follows. In order to make it 
self-contained, in Sec. [TT] we provide a brief summary of 
the RFA method for the structural properties of a multi- 
component HS mixture and state the result for the planar 
hard wall limit in which the concentration of one of the 
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species goes to zero while its diameter goes to infinity. 
The expUcit derivation of such a Umit is made in the 
Appendix. Section IIIII is devoted to the outhne of a re- 
cent version of DFT derived by one of us pi*]. This is 
followed in Sec. llVl bv a description of the simulation de- 
tails. In Sec. |V] we present the results derived with the 
two theoretical approaches as well as a comparison with 
the simulation data. The paper is closed in Sec. I VII with 
further discussion and some concluding remarks. 



II. THE RATIONAL FUNCTION 
APPROXIMATION METHOD 

We start by considering an {N + l)-component fluid 
of (additive) HS of diameters tT;, mole fractions Xi (with 
i = 0, 1, ... , N), and total bulk number density J). The 
partial bulk number density of species i is Pj = Xip and 
the packing fraction of the mixture is 77 = (7r/6)p(cr'^), 
where 



N 



(2.1) 



denotes the nth moment of the size distribution. 

According to the RFA for HS mixtures [H ill, the 
Laplace transform 



(2.2) 



(s) = / dre 'Wg.j (r) 



of rgij (r) , where gij (r) is the radial distribution function 
for the pair ij, is explicitly given by 



G.,(.s) = ^[L(s).B-i(. 



(2.3) 



where cTjj ~ (o-j +(Tj)/2 and L(s) and B(s) are {N + 1) x 
{N + 1) matrices given by 



(2)„2 



Bij{s) = (1 + as)Sij - Aij(s,, 



(2) 

+ipo{(Tis)aiLlj' 

In Eqs. 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



\m=0 



N 



Lf^ = A + AV, + 2A'a - A ^ Pfc^fcig , (2 



i-j' = Acry-t-^A'criCrj + (A-|-A'(Tj)a-iA(Ti^pj,(TfcL^^.\ 

(2.9) 



fc=0 



where 



A = 



27r 



A' = 



(2.10) 



In principle, the matrix L'^^ and the scalar a can be 
chosen arbitrarily without violating any basic condition 
[20l . [23 |. In particular, the choice L-^' = a = gives 
the PY solution [H, [2I]. One can go beyond this ap- 
proximation by prescribing given contact values gij(oij) 
and the associated thermodynamically consistent (re- 
ducedHsothermal compressibility x, what fixes L^^^ and 
a 0,12. Specifically, 



(2.11) 



while a is the smallest real root of an algebraic equation 
of degree 2(A^ + 1). Here, following the method of Ref. 
[131, we will take for gijicrij) the following extension to 
mixtures of the Carnahan-Starling-Kolafa contact value 
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37? (o-2)cr,crj 772(5-27/ + 2772) 



1-^7 2(1-77)' ('^') 



12(1-77)3 



772(1 + 77) / (CT2)CT,cr 



6(1-77)2 V(a3) a. 



(2.12) 



These contact values are thermodynamically consistent 
with Boublfk's equation of state for HS mixtures [2^ . 
The associated (reduced) isothermal compressibility is 



1 



677 (a)(a2} 

(1-77)2^(1-77)3 (a3) 
,27- 877 - 8772-1-4773 (0-2)3 



-'7 



3(1-77)^ 



(a3)2 



(2.13) 



Now we assume that the mole fraction of the spheres of 
species 7 = vanishes (xq — *■ 0) and that their diameter is 
infinitely larger than those of the other species (cro/di — > 
cx),7 > 1), in such a way that x^a^ <^C (cr") for 77 < 3. 
In that case, species does not contribute to either the 
total packing fraction 77 or the average values (cr), {a'^), 
and (cr^), i.e. 



(a") 



N 

1=1 



77 < 3. 



(2.14) 



fe=0 



Under those conditions a particle of species i = is seen 
as a planar hard wall by the A'^-component mixture made 
of particles of species 7 = 1, 2, . . . , A. By carefully taking 
the limits xq ^ and ctq — > 00 (with the constraint 
^o'^o (o'") for n ^ 3), it is proven in the Appendix 
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that the local density of particles of species i at a distance 
z from the wall is 



(2.15) 



where the function 7^(2) is the inverse Laplace transform 
of the function defined by Eqs. ((XT2)) . (f04)l - ((XT8)) . As 
said in connection with Eqs. ([^ - (j^ . (j^ . and 

p. lip , the PY theory for ji{z) is reobtained by setting 
a = 0. Our extended RFA theory is obtained by impos- 
ing Eqs. (|2.12p and (|2.13p and finding the corresponding 
value of a. Note that, within the RFA method, expres- 
sions for gij and x other than those of Eqs. (|2.12p and 
(|2.13p could equally be used. The choice we have made 
relies on the fact that these constitute well tested and 
reliable approximations. 



III. DENSITY FUNCTIONAL THEORY 

DFT is based on the property that, for a given inter- 
atomic potential, the grand potential fi, or, equivalently, 
the intrinsic Helmholtz free energy _F, is a unique func- 
tional of the one-body density profile pi{r) [l^- For an 
iV-component fluid mixture at given temperature T, to- 
tal volume V, chemical potential /i^, and external poten- 
tial V°^^{r) for each component, the equilibrium density 
profile Pi(r) minimizes the grand potential functional 

A' 

^ [{P^}] = F [{p,}] + J2 / drp,{r)[Vr\r) - f,,]. (3.1) 
i=i •' 

The ideal contribution to the intrinsic free energy func- 
tional is known exactly: 

N 

Fid [{p^}] = kBTj2 / drp.(r) {in [p.(r)A3] - 1} , 
i=i 

(3.2) 

where fc^ is the Boltzmann constant and A^ is the ther- 
mal de Broglie wavelength of species i. Thus, in the task 
of finding the appropriate free energy functional one can 
focus only on the excess part 



F,. [{p^}]=F[{p,}]-F,^ [{p,}] 



(3.3) 



Once the expression for the intrinsic Helmholtz free en- 
ergy functional is known, the density distribution is given 
explicitly by 

p.(r) = A-3 exp[cf )(r) + /3m, - pvr\r)], (3.4) 

where /3 = l/{kBT) and c\^\'r) is the one-body direct 
correlation function, 



SF.^ [{P,}] 
6pi{r) 



(3.5) 



Because the knowledge of the free energy functional 
provides full description of the model under study, prac- 
tically any implementation of DFT requires some explicit 



approximation for the functional Fc^. In the fundamen- 
tal measure theory [28j . the functional i^cx of a mixture 
of HS is assumed to take the form 

fox [{a}] = keT y"dr$({n,(r)}), (3.6) 

where the excess free energy density, $({nQ,(r)}), de- 
pends only on the system averaged fundamental measures 
of the particles 

N 

na(r) =Y.J drV,(r')^f ^(r - r'). (3.7) 



i=l 



The weight functions cjj-"^(r) characterize the geometry 
of particles. The minimal space of the weight functions 
is generated by the basis 



(r)=e(a,/2-r). 



(r)=5(r-a,/2). 



u;f)(r) = -5(r-a,/2) 
r 



(3.8) 



(3.9) 



(3.10) 



where 0(r) is the Heaviside step function and S{r) is 
the Dirac delta function. The other weight functions are 
proportional to those given by Eqs. (|3.8p - (|3.10p . namely 
^«(r) = (r)/27ra„ ^(r) = 4'\r)/naf, a;«(r) ^ 
wf^(r)/27ra,. 

In this work, we apply the following form of the free 
energy density [2l|: 



$ = 



, , , nin2 — rii • n2 / „ \ 
-no ln(l - 113) H h (71-2 - 3n3n2 • 112) 

8(1 - 713)2 ln(l - na) + 8723 - (15/2)7i^ + 2nl 



1087rni(l - 713)2 



(3.11) 



In the limit of a bulk fluid, the vectorial weighted den- 
sities rii and 112 vanish, tiq — 'p, ni = {'p/2){a), n2 — 
■7T{'p/2){a^), 7i3 = ?7, and therefore excess free energy den- 
sity becomes equivalent to that following from Boublfk's 
equation of state [2^ . In the low density limit the func- 
tional given by Eqs. (|3.6p and (|3.1ip is equivalent to that 
originally proposed by Rosenfeld, which underlies the PY 
compressibility solution. It is well known, however, that 
the PY compressibility equation of state overestimates 
the pressure at high densities. The functional presented 
here includes corrections which more accurately extrap- 
olate from low to high density states pl| . 

In the problem at hand, the external potential repre- 
sents the interaction with a planar hard wall, so that 



K'"' (r) - ^ for z > cr,/2, 
where z is the coordinate normal to the wall. 



/cxt/„^ _ J 00 for 2: < 0-^/2 



(3.12) 
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IV. SIMULATION DETAILS 

A binary mixture of HS confined in a box of dimen- 
sions {Lx, Ly, Lz) was simulated using the NVT Monte 
Carlo method. The particles were sorted out to cells and 
the linked-list method was used. Periodic boundary con- 
ditions were applied in the x and y directions, whereas 
flat hard walls were located at z = and z = L^. We 
chose L,j. — Ly = 10 and — 30 in units of ai (the 
diameter of the small spheres) for all the simulations. 

The initial parameters of the simulations correspond 
to the average densities and the average concentrations. 
The connection between average and bulk values is not 
initially known so that the bulk values are evaluated in 
the region of the simulation box z Lz/2, where no in- 
fluence of the walls is expected. It is worth mentioning 
that the bulk values can be alternatively determined us- 
ing two separate boxes connected by the same chemical 
potentials (one with and the other without hard walls) 
or using the Gibbs ensemble. However, we have found 
the simple NVT method (which does not involve parti- 
cle insertions) to be more efficient and less problematic, 
especially for asymmetric mixtures and/or for high den- 
sities. 

There are two common procedures as to how an ini- 
tial configuration can be prepared. The random shooting 
method starting with the species of the largest diameter 
fails to converge at high packing fractions. On the other 
hand, the box can be initially filled up with particles 
placed at the crystal configuration, usually an fee lattice. 
The disadvantage of the latter is the need of "relaxing" 
the system, a process that may take too long and may 
even turn to be untractable. We present a novel method 
which suffers none of the above mentioned problems and 
can be easily generalized to the model of an arbitrary 
number of components. First of all, both "large" and 
"small' particles are inserted randomly into the box with 
the only restriction of no overlapping of the like parti- 
cles. Then the system is left to propagate keeping the 
interaction between the like particles hard, whereas the 
penetrability of the unlike particles is restricted by the 
potential 

"^J^'^^-jo for rXTy, 

where the value a 1 of the constant turned up to be 
convenient. Choosing the initial reduced temperature 
(typically T* = ksT/e — 0.05), the system was gradually 
cooled down. We found this algorithm to be extremely 
efficient in creating a non-overlapping configuration (typ- 
ically within one minute using a standard PC). 

The system is subsequently equilibrated over more 
than 10^ MC steps and the bulk densities and density 
profiles are averaged over another 4 x 10^ moves. The 
average runs were divided into 20 subaverages to esti- 
mate the standard deviations. 



V. RESULTS 

In this Section we report the results we have obtained 
with the previous approaches, namely the RFA method, 
DFT, and simulation, for the density profiles of binary 
HS mixtures in the presence of a planar hard wall. In 
order to test the theories under extreme conditions, we 
consider cases in which the sizes of the two components 
in the mixture are disparate but both occupy a similar 
volume. 

For the sake of illustration, in Figs. [lH3] we present 
the density profiles for both components of nine binary 
mixtures with a common size ratio (72 /ci = 3 and a total 
(bulk) packing fraction r] ~ 0.2 (Fig. [1]), -q ~ 0.3 (Fig. 
[5]), and Tj ~ 0.4 (Fig. [3]). For each packing fraction, three 
different compositions have been considered: X2 — 0.06, 
so that 772/?7i — 1.7 (top panels); X2 — 0.04, so that 
?72/?7i — 1.1 (middle panels); and X2 — 0.02, so that 
ri2/rii ~ 0.55 (bottom panels). Here, rji = (7r/6)pjCrf {i = 
1, 2) denotes the partial (bulk) packing fraction of species 
i. Since rj and X2 are measured in the bulk, they present 
minor deviations with respect to the imposed average 
values in each case. Apart from the simulation data and 
our theoretical approaches, the figures also include the 
PY results. 

Without any doubt, the DFT is the one that produces 
the best overall agreement. As far as the RFA method is 
concerned, one finds that it also does a very good job in 
general, being particularly accurate not only at the con- 
tact value (which is of course an input in this approach), 
but also caters very well for the second maximum and 
the rest of the oscillations. These nice features worsen 
in the vicinity of the first minimum, particularly for the 
density profile of the bigger species at 7y ~ 0.4, where 
the RFA method may lead to the unphysical prediction 
of a negative value. As expected, the PY theory yields 
poorer contact values and gets worse as the total pack- 
ing fraction increases, yielding a negative value for the 
first minimum of P2{z) at 77 ~ 0.4 and X2 — 0.02. In any 
case, it indeed accounts for the oscillatory character of 
the profiles. 

Figures [THSl show that a rich and complex structure 
of the local densities pi{z) and P2{z) appears for all the 
cases considered. This is related to the fact that rji ^ rj2 
and so both species compete for the available volume. As 
the total density increases so does the nonuniformity of 
the density profiles, as reflected by the contrast between 
the values at contact and at the first minimum, especially 
in the case of the large spheres. Moreover, as the mole 
fraction X2 of the large spheres decreases, the character- 
istic wavelength of the oscillations decreases. To better 
understand this phenomenon, imagine a mixture with a 
mole fraction X2 such that 771 <C 772 • In that case, the 
large spheres (species 2) are practically unaffected by the 
presence of the small ones, so the local density P2{z) is 
the same as that of a single component fluid at the same 
packing fraction and oscillates with a characteristic wave- 
length of the order of 0-2 ■ As for the small spheres (species 
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FIG. 1: (Color online) Density profiles close to the wall for both species in the binary HS mixtures with cr2/o"i = 3 having 
a total packing fraction close to ?7 = 0.2. Top panels: rj = 0.208 and X2 = 0.063; middle panels: i] — 0.207 and X2 ~ 0.042; 
bottom panels: rj — 0.206 and X2 = 0.021. Solid lines: DFT; dashed lines: RFA results, dash-dotted lines: PY theory. The 
symbols are the results of simulation. The insets provide a wider range and more details of the resulting structure. 



1), their density profile is dominated by the presence of 
the large spheres and so the wavelength of pi{z) is also of 
the order of a2- In the other extreme situation, namely 
when ■i]2 Tji, we have the opposite situation: the small 
spheres behave as a single component fluid and enslave 
the density profile of the large spheres, so that both pi{z) 
and p2iz) oscillate with a characteristic wavelength of the 
order of ai. The interplay between both extreme cases 
occurs when rji ~ 772, giving rise to the superposition of 



both length scales, as observed in Figs.[TH3l To illustrate 
the transition from 771 ^ 772 to 772 ^ in Fig.|3]we plot 
the density profiles predicted by the RFA for mixtures 
with (T2/cri — 3, T] — 0.3, and 2:2 = ^ — 0.27 (which 
corresponds to 771/7/2 — 0.1), X2 = ^ — 0.036 (which cor- 
responds to rii/ri2 = 1), and X2 = — 0.0037 (which 
corresponds to T]i/rj2 — 10). A ratio 1:10 in the par- 
tial packing fractions is enough to make the species with 
the largest packing fraction to behave almost as a single 
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FIG. 2; (Color online) The same as Fig. \T\ but with a total packing fraction close to ?7 = 0.3. Top panels: rj = 0.309 and 
X2 = 0.062; middle panels: rj = 0.308 and X2 ~ 0.037; bottom panels: rj = 0.307 and X2 = 0.020. 



component fluid. In fact, Fig. [5] shows that the curve 
of P2{z)/p2 in the case X2 = — 0.27 (?7i/?72 = 0.1) 
and the curve of pi{z)/'pi in the case X2 = 571 — 0.0037 
(?7i/'72 = 10) look Hke very similar when each one is plot- 
ted as a function of the corresponding scaled distance 
z/ai. It must be noticed that, as observed in Fig.O the 
case 771/772 — 10 is actually closer to the single component 
fluid than the case 771/772 = 0.1. In the latter case one 
has X1/X2 = 2.7, Xiai/x2cr2 — 0.9, and Xiaf/x2cr2 = 0.3, 
so that the influence of the species 1 is not entirely negli- 
gible. On the other hand, the ratios are X2/X1 = 0.0037, 



X2cr2/xiai = 0.011, and a;2cr2/^i'^i — 0.033 in the case 
m/v^ = 10. 

VI. CONCLUDING REMARKS 

The results presented in the preceding Section deserve 
some further discussion. To begin with, inspired by the 
work of Noworyta et al. 0] , we have revisited the prob- 
lem of determining the structure of binary HS mixtures 
in the presence of a planar hard wall using three of the 
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FIG. 3: (Color online) The same as Fig. [T] but with a total packing fraction close to 77 = 0.4. Top panels: X2 = 0.062 and 
T] = 0.401; middle panels: X2 = 0.036 and rj = 0.406; bottom panels: X2 = 0.019 and 77 = 0.404. 



different approaches that have been proposed in the ht- 
erature to deal with it. Concerning the RFA method, its 
main advantage is that of allowing for a completely ana- 
lytical description, as it also occurs with the PY theory, 
avoiding at the same time the thermodynamic inconsis- 
tency problem present in the latter. It is fair to say that, 
under the rather extreme conditions that we tested it, 
this method is a reasonable compromise between accu- 
racy and simplicity. Whether the difficulties associated 
with the prediction of a negative first minimmn of p2{z) 



at 77 ~ 0.4 and X2 S 0-04 could be avoided by the intro- 
duction of additional parameters in the method remains 
to be assessed. It should be mentioned in passing that 
for mixtures with 77 ~ 0.3 this unphysical behavior is not 
predicted, even for X2 —>■ 0. Moreover, in the case of pure 
ternary mixtures, the agreement between the results for 
the structural properties of the RFA method and simu- 
lation has been shown to be highly satisfactory [2^. Fi- 
nally, although wc have illustrated the results for the case 
of binary mixtures, a further asset of the development we 
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FIG. 5: (Color online) Density profiles close to the wall in the 
binary HS mixtures with a^joY = 3 having a total packing 
fraction r\ = 0.3, as predicted by the RFA. Solid line: pi (£) / pi 
versus z/a\ for X2 — ^ 0.0037 {rj\/ri2 = 10); dashed line: 
P2{z)/p2 versus z/a2 for X2 = ^ ~ 0.27 {ni/v^ = 0.1). The 
dotted line corresponds to the single component case. 



have presented applies in principle to any multicompo- 
nent mixture of HS near a planar hard wall. On the 
other hand, the excellent performance of the DFT in the 
approximation introduced in Ref. plj has been already 
pointed out. However the merit of the proposed excess 
free energy functional cannot be overlooked and might 
be useful for other purposes as well. As a final point, it 
should be stressed that the simulation method reported 
in this paper is also novel and allowed us to obtain re- 
sults for situations that were problematic ten years ago. 



In particular, it allowed us to confirm the assertion made 
by Roth and Dietrich concerning the non existence 
of the anomaly reported by Noworyta et al. [8|] for the 
mixture with size ratio 1:3. Our hope and expectation 
is that this method proves useful as well for other inter- 
esting systems, some of which we plan to examine in the 
future. 



Before closing this paper and for the sake of setting a 
wider perspective for the results we have presented, a few 
further comments are in order. A simple and yet realistic 
model of colloidal dispersions consists of a highly asym- 
metric binary hard-sphere mixture in which the large 
spheres stand for the colloidal particles and the small 
spheres represent solvent or polymer molecules. There- 
fore, our results may also be applied in these systems, 
which are nowadays easily amenable for experimental ex- 
amination, and open up the possibility of investigating 
solvation forces and other interesting physical phenom- 
ena associated with particles of different sizes competing 
for interfacial positions in the presence of walls. We plan 
to pursue some of these issues in the future. On the other 
hand, the possible extension of our work to deal with non- 
hard particles remains to be explored. In the case of the 
RFA method, structural properties of sticky hard-sphere 
mixtures and other systems have already been derived 
[2^ . As for the DFT approach, modifications would be 
required. Thus it seems that much work remains to be 
done before such extensions become a reality. Again, 
we plan to explore these and other possibilities in future 
studies. 
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APPENDIX: THE WALL LIMIT IN THE RFA 

The limit — > implies that the row i = of the 
(iV-l- 1) X (A^-l- 1) matrix A defined by Eq. ((2?6| vanishes, 
so that 

Ay = ly (1 - ,5,o)(l - Sjo) + Aio{l - SM)SjQ. (A.l) 

Here, Aij is a projected N x N matrix with i,j>l. In 
fact, Aij is the matrix corresponding to an iV-component 
mixture in the absence of species 0. Inserting (jA.l|l into 
Eq. (|2.5p . we have 

Bij = Bij{l-Sto){l-Sjo) + {l+as)StoSjo-Aio{l-5io)Sjo, 

(A.2) 

where Bij is the N x N matrix 

B,j = (1 + as)6ij - Aij , ij>l. (A.3) 

It can be checked that the {N -I- 1) x {N + 1) inverse 
matrix B^^ is 

(B-^) = (b-i) (1 - d,o){l - S,o) + (1 + as)-'6,o5,o 



-Ci{l - 5io)5jo, 



(A.4) 



where B ^ is the N x N inverse matrix of B and the 
elements C, are 



1 ^ 



(A.5) 



Insertion of Eq. (|A.4|) into Eq. (|2?3)) yields 



2tts- 



N 

^e---Gy = ^L,fc(B-i) (l-5,o) + T 



Lio 



k=l 



N 



LikCkSjf). 



(A.6) 



In particular, if i,j > 1, 



N 



fc=l 



Ik] 



hi > 1, 



(A.7) 



if > 1- Equation (|A.6p implies that, as 



where we have taken into account that | B 
(B- _ 

expected, the iV-component mixture is not affected by 
the presence of the species in the infinite dilution limit 
xq 0. On the other hand, setting i > 1 and j = in 
Eq. (IA.6p . we have 



Lio 



1 



2Tr^ 
l + a 



N 



Ve'^--G„Ao, i>l, 



(A., 



where use has been made of Eqs. (jA.5p and (|A.7p . The 
cross function Gio (with i = 1, . . . , N) is related to the 
spatial correlation between the diluted species and the 
species i of the true A'^-component mixture. We see from 
Eq. (|A.8[) that Gio is expressed in terms of the matrix 
Gij of the A^-component mixture and the cross elements 
Lio and Ajq. 

Let us now introduce the shifted radial distribution 
function 



li{z) = gmiz + (To/2), 



(A.9) 



where z > ai/2 represents the distance from the center of 
a sphere of species i to the surface of a sphere of species 
0. In Laplace space. 



G,o(s) 



-CTos/2 



CTq 



r.(s)-rK.) 



where 



r.(s) 



dze '*''7i(z) 



(A.IO) 



(A.ll) 



ai/2 



is the Laplace transform of "/i{z) and r,^(s) = dri(s)/ds. 

Thus far, the diameter (Tq is arbitrary as long as Eq. 
(I2.14P is satisfied. Now we take the wall limit (do oo). 
In that case, the function 'jiiz) has a clear meaning as the 
ratio between the local density of particles of species i at 
a distance z from the wall, Pi{z), and the corresponding 
density in the bulk, pj = pi(oo). In the wall limit r^(s) 
can be neglected versus ao Ti{s) inEq. (jAlOl) . Therefore, 

T,{s) - 2 lim ao-'e''«'/^G.,Qis) 



-CTiS/2 



Ao(s) 



N 

xY^e^ 



7rs2(l -|- as) 
""Gy(s)Ajo(s) 



as 



i > 1,(A.12) 



k=l 



where in the last step use has been made of Eq. (|A.8|) 
and 

Ao(s) = lim cro"^Ljo(s), Ao(s) = lim CTq ^Ajo(s). 

CTq — >00 (To — 

(A.13) 
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From Eqs. ([23), jM]), dH]), (HH), and (|2lT|) one gets 



AT 



4? = A' - TraA ^ Pjajgjo{ajo), (A.15) 



A A' A 



(A.17) 



A^o{s) Pi (p2(cr^s)cr,^£|°^ +¥5l(crjs)cr,^/:,|o^ 



(2) 



(A.18) 



It must be noted that the parameter a appearing m Eqs. 
(IA.12P and (|A.15p - (|A.17p is independently obtained from 



the RFA solution for the iV-component mixture. Regard- 
ing the contact values (?io(cio)i they are obtained by tak- 
ing the limit cto ^ cxd in Eq. p.l2p . namely 



'JiO 



(^'•o) = Z + 



3ry (cr2) r/2(5-2?7 + 2r/2) 



1-,, (1-^f (a3} ^ 3(1-7^)3 



-cr-i 



(A.19) 



Taking into account that G.y (s) = s ^ + 0(s") for 
small s, it is possible to prove that 



lim sVi(s) — 1, 



which implies the physical condition 
lim 7i(z) = 1. 



(A.20) 



(A.21) 
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